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ABSTRACT 

We present deep optical and X-ray follow-up observations of the bright unassociated Fermi-LAT gamma- 
ray source 1FGL J131 1.7-3429. The source was already known as an unidentified EGRET source (3EG J1314- 
3431, EGR J1314-3417), hence its nature has remained uncertain for the past two decades. For the putative 
counterpart, we detected a quasi-sinusoidal optical modulation of Am ~ 2 mag with a period of ~1.5 hr in 
the Rc, r' and g' bands. Moreover, we found that the amplitude of the modulation and peak intensity changed 
by > 1 mag and ^0.5 mag respectively, over our total six nights of observations from 2012 March and May. 
Combined with Swift UVOT data, the optical-UV spectrum is consistent with a blackbody temperature, kT 
~ 1 eV, and the emission volume radius ~ 1.5 x 10 4 dkpc km (dkpc is the distance to the source in units 
of 1 kpc). In contrast, deep Suzaku observations conducted in 2009 and 2011 revealed strong X-ray flares 
with a lightcurve characterized with a power spectrum density of P(f) oc /-2.o±o.4 but the folded X-ray light 
curves suggest an orbital modulation also in X-rays. Together with the non-detection of a radio counterpart, 
and significant curved spectrum and non-detection of variability in gamma-rays, the source may be the second 
"radio-quiet" gamma-ray emitting milli-second pulsar candidate after 1FGLJ2339. 7-0531, although the origin 
of flaring X-ray and optical variability remains an open question. 
Subject headings: gamma rays: stars — pulsars: general — X-rays: general 



1. INTRODUCTION 



The Large Area Telescope (LAT; lAtwood et al. 2009]) on- 
board the Fermi Gamma-ray Space Telescope is a succes- 
sor t o EGRET onboard th e Compton Gamma-ray Observa- 
tory (lHartman et aljfl999l) . with much improved sensitivity, 
resolution, and energy range. The second Fermi-LAT cat- 
alog, b ased on the first 24 months of all-sky survey data 
(2FGL; lNolan et aLl l2012). provides source location, flux and 
spectral information, as well as light curves on month time 
bins for 1873 7-ray sources detected and characterized in 
the 100 MeV to 100 GeV range. Thanks to their small lo- 
calization error circles (or ellipses) with typical 95% confi- 
dence radii, ~ 0.1° — 0.2°, for relatively bright sources, 
69 % of the 2FGL sources are reliably associated or firmly 
identified with counterparts of known or likely 7-ray pro- 
ducing sources. In particular, more than 1,000 sources are 
proposed to be associated with AGN (of mainly the blazar 
class) and 87 sources with pulsars (PSRs), including 21 mil- 
lisecond pulsars (MSPs) which are a n ew category of 7- 
ray sources discov ered with Fermi-LAT dNolan et al J 120121 : 
lAbdo etal.ll2009ah . Other sources, albeit of a relative mi- 
nority compared to AGNs and PSRs, also constitute im- 
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portant categ ories of new GeV s ources like supernova rem- 
nants rSNRs: lAbdolt^l2009b[ l2010albTch . low-mass/high- 
mass binari es ( Abdo et alj l2009ctfdlleh. pulsar win d nebula 
dAbdo et alj|2010dlleh . one nova dAbdo et al 1120101 , normal 
and starburst galaxies dAbdo et alJ l20103lhl) . and the giant 
lobes of a radio galaxy (Abdo et al.l l2010il) . 

Despite such great advances in the identification of Fermi- 
LAT sources, 575 (31%) sources in the 2FGL catalog still 
remain unassociated. Note that a substantial fraction of the 
unassociated sources (51%) have at least one analysis flag 
due to various issues, w hile only 14% of the associated 
sources have been flagged dNolan et alJ l2012). This may sug- 
gest some of unassociated sources are spurious due to com- 
plexity/difficulty of being situated in a crowded region near 
the Galactic plane. Nevertheless, many of them are bright 
enough to be listed in the one year Fermi-LAT catalog ( 1 FGL; 
I Abdo et alj|2010ih and some of them are even listed in the 
bright source lis t based on the first 3 months of data (0FGL; 
Abdo et al. 2009f). By comparing the distribution of associ- 
ated and unassociated sources in the s ky, a number of in terest- 
ing features in the map were reported dNolan et al.ll2012l) . For 
example, (1) the number of unassociated sources decreases 
with increasing Galactic latitude, (2) the number of unasso- 
ciated sources increases sharply below Galactic latitudes, |6| 
< 10°, and (3) the fraction of sources with curved gamma- 
ray spectra among the unassociated sources is greater (28%) 
than the fraction of curved spectra sources among the asso- 
ciated sources (16%). Further extensive studies based on a 
statistical approach in an effort to correlate their gamma-ray 
properties with the AGN and P SR populations was pre sented 
for 1FGL unassociated sources ( Acker mann et alj l2012). 

In this context, 1FGL J131 1.7-3429 (or 2FGL name, 
2FGL J131 1.7-3429) is a classical unassociated gamma-ray 
source situated at high Galactic latitude (I = 307.6859°, b 
= 28.1951°), and was first discovered by EGRET about 20 
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years ago as 3EG J1314-3431 dHartman et al.lll999h or EGR 
J1314-3417 dCasandiian & Grenieri 120081) . The source was 
also reported by ferm/-LAT in the OFGL list with a gamma- 
ray flux F .i_20GcV = (11.7 ±1.1) x 10~ 8 photons cm~ 2 s _1 , 
which is marginally consistent with the gamma-ray flux deter- 
mined by EGRET, i*o.i-20GcV = (18.7 ±3.1) x 10~ 8 photons 
cm~ 2 s _1 , within the 2a level. In the 2FGL catalog, the de- 
tected significance of 1FGL J131 1.7-3429 is 43.1 a, which 
is one of the brightest sources with an UNID flag. Based on 
the one-month binned > 100 MeV gamma-ray light curve of 
1 FGL J 1 3 1 1 .7-3429 over two years of data as published in the 
2FGL catalog, no statistically significa nt variability was ob - 
served ( Variabilityjndex = 19.09; iNolan et all (120121) ). 
The gamma-ray spectrum is significantly curved with curved 
significance SlGNlF_CURVE of 6.33 ( INolan et al.ll2012h .PI 

The first X-ray follow-up observation of 1FGL J1311.7- 
3429 was conducted as a part of Suzaku X-ray observa- 
tions of 11 unidentified Fermi-LAT object s at high Galac- 
tic la titude, |6| > 10° dMaeda et alJ 1201 U iTakahashi et al.l 
120121) . The X-ray source associated with 1FGL J131 1.7-3429 
showed a very rapid X-ray flare with the count rate chang- 
ing by a factor of 10. Subsequent Chandra ACIS-I (2010 
March 21 for a 19.87 ks exposure, obsID 11790) and Swift 
XRT observations (2009 February 27 for a 3.34 ks exposure, 
obsID 31358; see also Table 1) confirmed that the brightest 
X-ray source within the Fermi-YKY error ellipse is the most 
credible counterpart and that the X-ray source is also variable 
on month-to-year timescales. The unabsorbed X-ray flux ob- 
served with Chandra was 1.03x 10 -13 erg cm -2 s _1 in the 
0.5—8 keV band, with a differential photon spectral index, Y 
= 1.26±g;lf (ICheung et alJl2012h . 

Motivated by the initial X-ray results, we conducted further 
deep observations of 1FGL J131 1.7-3429 with Suzaku, to- 
gether with deep optical observations using a 105 cm Ritchey- 
Chretien telescope (<?', Rc, and Ic bands) at the Ishigakijima 
Astronomical Observatory (IAO) in Japan, as well as the lm 
telescope (r' band) at Lulin Observatory in Taiwan. In Sec- 
tion 2, we describe the details of the Suzaku observations and 
optical observa t ions and data reduction procedures. Very re- 
cently, [Romani| (|20T2|) reported quasi-sinusoidal optical mod- 
ulation of this source with a 1.56 hr (5626 s) period, suggest- 
ing that the source is another black- widow-type M SP like that 
recent ly discovere d for 1FGL J2339.7-0531 dRomani & Shawl 
l2^lKongetani20l2l) . Our paper confirms some of those 
optical findings for 1FGL J131 1.7-3429, plus provides re- 
sults from multiple epoch optical monitoring between 2012 
March and May and completely new X-ray data based on a 
long Suzaku observation conducted in 2011 together with 
our previously published archival 2009 data. The results of 
these observations are given in Section 3. Based on our new 
observational data in optical and X-ray, and various observed 
gamma-ray parameters compiled in the 2FGL source catalog, 
we support that 1FGL J131 1.7-3429 could be a "radio-quiet" 
gamma-ray emitting MSP candidate like 1FGL J2339.7-053 1 . 

8 http://heasarc.gsfc.nasa.gov/FTP/fermi/data/lat/catalogs/source/lightcurves/ 
2FGL J13 1 Id7m34 29 Jcpng 

9 As detailed in Nolan et al. 1 2012), the VARIABILITYJNDEX is an in- 
dicator of the flux being constant across the full 2-year period. A value of 
VARIABILITYJNDEX > 41.6 is used to identify variable sources at a 99% 
confidence level. Similarly, the SlGNIF_CURVE parameter is an indicator of 
the spectrum being curved, by comparing the likelihood values calculated for 
a LogParabola and a single power-law function. SIGNIF_CURVE is distributed 
as x 2 with one degree of freedom, and SlGNIF.CURVE > 16 corresponds to 
4 a significance of curvature. 
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FIG. 1 . — Tri-color optical image of 1FGL J131 1.7-3429 constructed from 
data taken in the g', Rc, and Zc-bands by the IAO 1 .05 m telescope (Table 1). 
The four reference stars used in the analysis of the IAO data, are denoted as 
Ref 1-4, while two refrence stars used in the analysis of the LOT data are 
shown as boxes. 

The variable optical/X-ray source is posited as the counterpart 
to the gamma-ray source and throughout, we refer to it simply 
as 1FGLJ131 1.7-3429. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Optical/UV 

As discussed in C heung et al.l (120121) . the brightest X- 
ray source within the Fermi-LAT error ellipse, CXOU 
J131 145.71-343030.5, is the most credible counterpart and 
exactly the one dete cted in our previ ous (A04; below) 
Suzaku observation dMaeda et al.1 1201 lb . The source was 
also detected with the Swift XRT (observation ID 31358). 
Within the Chandra error circle (0.6" at 90% level), there 
is a R = 17.9 mag, B = 20 .5 mag star in the USNO A-2.0 
catalog dMonet et alj|2003al) . that is 0.56" apart from CXOU 
J131 145.71-343030.5. The same star is also listed as R = 
18.8 mag, B = 21.0 mag in USNO-B1.0 catalog, suggesting 
a hint of temporal variability (a typical uncertai nty of these 
measurements is ~ 0.3 mag ; Monet et al. 2003a). The same 
optical source is also seen in the DSS (Digital Sky Survey), 
but its optical magnitude is unclear. Moreover, Swift UVOT 
observations taken simultaneously with the XRT detected the 
source. Our analysis using archival Swift UVOT data in- 
dicate: v > 19.51 (upper limit only), b = 20.10±0.30, u = 
20.77±0.27, uvwl = 21.70±0.31, uvm2 = 21.58±0.24, and 
uvw2 = 22.05±0.22. 

We made further deep follow-up observations of the field 
of CXOU J131 145.71-343030.5, centered at (RA, DEC) 
= (197.940400°, -34.508306°), with the 105 cm Ritchey- 
Chretien telescope at the Ishigakijima Astronomical Obser- 
vatory in Japan. These observations were obtained on 2012 
May 25 and started at 12:31:11:13 and ended 15:41:11:63 
(UT). The telescope in IAO is equipped with a tricolor cam- 
era that performs simultaneous imaging in the SDSS-g' (here- 
after, g'), Rc and Ic bands. The total net exposure amounts 
to 8400 s (300 s x 28 frames; see Table 1 for the observation 
log). All images were flat field and bias corrected. The ab- 
solute magnitudes were calibrated against the four reference 
stars shown in Figure 1. In the Rc and Ic bands, the ref- 
erence star magnitudes were based on the NOMAD catalog 
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TABLE 1 

Observation log of 1FGL J131 1.7-3429 analyzed in this paper 



Obs Start (UT) 


Obs End (UT) 


Observatory 


Band 


Exposure 


2009-08-04 04:56:35 


2009-08-05:07:18:14 


Suzaku XIS 


X-ray 


33.0 ks 


2011-08-01 16:48:20 


2011-08-03:17:40:15 


Suzaku XIS 


X-ray 


65.2 ks 


2009-02-27 18:41:59 


2009-02-27:22:04:57 


Swift UVOT 


UV 


276/276/276/552/742/1 106 a 


2012-03-24 17:45:37 


2012-03-24 19:38:35 


LOT lm 


r' 


5 min X 21 


2012-03-25 17:48:26 


2012-03-25 19:26:33 


LOT lm 


r' 


5 min X 12 


2012-03-26 17:25:14 


2012-03-26 18:29:59 


LOT lm 


r' 


5 min X 12 


2012-04-12 16:35:08 


2012-04-12 18:24:30 


LOT lm 


r' 


5 min X 21 


2012-05-24 13:00:49 


2012-05-24 15:34:44 


LOT lm 


r' 


5 min X 29 


2012-05-25 12:31:11 


2012-05-25 15:41:11 


IAO 105cm 


g', Rc, Ic 


20 min X 7 



NOTE. — a : Swift UVOT exposures for the v/b/u/uvwl/uvm2/uvw2 bands in seconds. 
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FIG. 2.— Multi-band optical light curves of 1FGL J131 1.7-3429 observed 
with the IAO 1.05 m telescope on May 25 (g', Rc and Ic bands). Data are 
binned by 1200 s (300 s X 4 frames) for the IAO data. The source is not 
detected in the Zc-band and 3 a upper limits are provided. 



dMonet et al.l l2Q03bb . Because the target area was not cov- 
ered by the SDSS, we employed a sy stem conver sion formula 
of USNO-B Be and Rc magnitudes (ISesarll2006l) for the four 
reference stars to obtain g' magnitudes. Moreover, the ob- 
served magnitudes were corrected for Galactic extinction us- 
ing A g , = 0.247 mag, A Rc = 0.167 mag, and A Ic = 0.120 
mag. 

We also conducted optical SDSS-r' (hereafter r') band 
monito ring observations w ith the Lulin One-meter Telescope 
(LOT; iHuang et al.l l2005). Observations were conducted on 
five nights between 2012 March and May (see Table 1). Pho- 
tometric images with 300 s exposure were obtained using the 
PI1300B CCD camera. We performed the dark-subtraction 
and fiat-fielding correction using the appropriate calibration 
data. For these LOT data, the four reference stars used in the 
IAO photometry were saturated in the detector so could not 
be used. Instead, the LOT photometric results are presented 
as differential magnitudes against two other fainter reference 
stars (boxed in Figure 1). 

2.2. Suzaku XIS 

As noted above, the first X-ray follow-up observation of 
1FGL J131 1.7-3429 was conducted in 2009 as a part of AO- 
4 Suzaku program (PI: J. Kataoka; OBS.ID 804018010) 
aimed at observing an initial four out of 11 unidentified 
Fermi-L AT objects a t high Galactic latitu de, \b\ > 10° (Maeda 
et al. 201 1 ; see also Takahashi et al.l2012h . To further investi- 
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FIG. 3. — (top) Temporal variations in the differential r'-band magnitudes 
derived by comparison with a reference star in the field, observed with the 
Lulin lm telescope from 2012 March to May. The r'-band magnitude of 
reference star remains constant within Am = 0.068 mag. Note that the mod- 
ulation profile differs largely among five nights of observations, (bottom) 
Folded light-curve of the differential r'- band mag nitudes with a best-fit pe- 
riod of 1.56278 hr (5626 s) proposed by Romani (2012). The phase zero is 
defined as MJD 56010.76808. Note that both the amplitude of the modula- 
tion and the peak intensity changed by > 1 mag and ~ 0.5 mag respectively, 
over the five nights of observations. 



gate the nature of the detected variable X-ray counterpart, we 
conducted the second Suzaku observation of 1FGL J131 1.7- 
3429 in 2011 as a part of AO-6 program (PI: J. Kataoka, 
OBSJD 706001010). The observation started at 2011 Aug 
01 16:48:20 and ended at Aug 03 17:40:15. The total ex- 
posure amounted to 65.2 ks, and is almost twice as long as 
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FIG. 4.— Suzaku XIS FI(XIS0+3) images of the 1FGL J1311.7-3429 re- 
gion in the 0.5—10 keV photon energy range using the data during the first 
20 ks (top; flare) and the last 74 ks (bottom; after the flare) in the 2009 ob- 
servation obtained in A04. The image shows the relative excess of smoothed 
photon counts (arbitrary units indicated in the color bar) and is displayed 
with linear scaling. The thick solid ellipse denotes the 95 % position error 
of 1FGL J131 1.7-3429 in the 2FGL catalog, while the dashed ellipse shows 
that reported in the 1FGL catalo g. Note the f aint source to the south presented 
(and denoted as src B) in Maeda et al. (201 1) is outside the 2FGL error circle. 

that obtained in AO-4 (see Table 1). For both the AO-4 and 
AO-6 data analysis, we excluded the data collected during the 
time and up to 60 s after the South Atlantic Anomaly (SAA), 
and excluded data corresponding to less than 5° of the an- 
gle between Earth's limb and the pointing direction. More- 
over, we excluded time windows during which the spacecraft 
was passing through the low cut-off rigidity (COR) of below 6 
GV. We set the same source region to within a 1 ' radii around 
the respective X-ray flux maximum and the selection criteria 
for the data analysis for the two datasets were completely the 
same. Although Suzaku also carries a hard X-ray detector 
(HXD), consisting of the PIN and GSO, hereafter, we do not 
use the data because the source is too faint to be detected with 
HXD/PIN or GSO. 

3. RESULTS 

3.1. Optical/UV 

Figure 1 shows a multicolor image of 1FGL J131 1.7-3429 
constructed from the IAO g' (blue), Rc (green), and Ic (red) 
data from 2012 May 24. The optical counterpart of the X- 
ray source is clearly detected in the image, but apparently the 
source is rather "blue" compared to reference stars and nearby 
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FIG. 5.— Multi-band Suzaku X-ray light curves of 1FGL J131 1.7-3429 
obtained during A04 (2009) and A06 (2011) observations; 0.4-1 keV, 1-2 
keV, 2-4 keV, and 4-8 keV (from top to bottom). The XIS-0, 1, 3 data are 
summed. 

galaxies. Aperture photometry yielded average magnitudes 
of the source, g' = 20.97 ± 0.13, Rc = 21.17 ± 0.16, and 
Ic < 18.38 (3 a upper limit). We also checked the temporal 
profile of the optical emission for the IAO data (Figure 2). 
Relative photometry for the IAO images against the four field 
stars shows large amplitude (Am ~ 2 mag) quasi-sinusoidal 
modulation in the q' and Rc band s with a a timescale of 1.5 
hr, as suggested by Romani (2012). 

To search for possible long-term variability, Figure 3 (top) 
compares the r'-band light curve reconstructed only from 
LOT data between 2012 March and May. The differential 
magnitudes are derived by comparison with two reference 
stars in the field (see, Figure 1). Note, the r'-band magni- 
tude of the reference stars remain constant within Am = 0.068 
mag over the five nights of observations. Interestingly, the 
amplitude of modulation and the lightcurve profile seems to 
have changed among the five nights of observations. Specif- 
ically, a clear modulation of Am ~ 2 mag is visible in the 
May 24 data, whilst Am ~ 1 mag in the March 24 data, and 
almost unseen (Am < 0.2 mag) in the March 26 data. More- 
over, the peak magnitude differs by Am ~ 0.5 mag among 
the five nights, much larger than the fluctuations in the mag- 
nitudes of reference stars. Figure 3 (bottom) shows the folded 
light-curves of differential r'-band magnitud es with a best-fi t 
period of 1.56278 hr (5626 s) proposed by iRomanl (120121) . 
Phase zero is defined here as MJD 56010.76808, so that the 
time of the observed minimum r'-band magnitude in the May 
24 data is set at orbital phase <fi = 1.0. The peak r'-band 
magnitude is around ~20.5 for five nights of data from 2012 
March and April. Note, this is exactly consistent with what 
has been observed with IAO on May 25 (Figure 2; right). 
To further investigate the temporal variability/flaring in the 
UV data, we also re-analyzed the archival Swift UVOT data. 
However, the UVOT exposures for each filter (b, u, uvwl, 
uvm2, uvw2) were too short to search for variability (see Ta- 
ble 1). 

3.2. Suzaku XIS 

During the AO-4 observation, significant X-ray variability 
was detected at the beginning of the ob servation, where th e 
count rate changing by a factor of 10 (Maeda et al. 2011). 
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FIG. 6. — (top) Energy dependence of X-ray variability of 1FGL J1311.7- 
3429. The variability parameter, excess variance was calculated for the total 
exposures in A04 (2009; red) and A06 (2011; blue) in four energy bands 
(see Figure 5). (bottom) Normalized PSD (NPSD) calculated from the X- 
ray light curves of 1FGL Jl 31 1.7-3429. The dotted line shows the best-fitting 
power-law function of oc/ -20 . 
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FIG. 7.— The folded X-ray light curves of 1FGL Jl 3 11.7-3429 recon- 
structed from the XIS data taken in Suzaku A04 (top; red) and A06 
(bottom) observations. All the XIS0, 1, 3 data are summed in the energy 
range of 0.4- 10 keV. The phase zero is defined as MJD 56010.76808. The 
vertical axis represents the source counts divided by the average source in- 
tensity in the frame. Dashed l ine represen ts the folded X-ray light curves 
of a nearby field source (src B in Maeda et al. 20 1 1 ) to confirm that observed 
excess is not due to artifacts caused by orbital gaps in the Suzaku data. 



TABLE 2 

Fitting Parameters of Suzaku data for the Power-law Model 



Parameter 


A04 
Value and Errors™ 


A06 
Value and Errors™ 


JV H (10 20 cm" 2 ) 

r 

Flux (0.5-2 keV) 6 
Flux (2-10 keV)° 


noo+O.lS 
"•■"-0.13 
, 7 , +0.18 

99+0-27 

, 7 o+0.16 


< 0.07 
1 26+ - 08 

n no+0.08 

2.96±0.15 


X A (dof) 
P(X 2 ) 


17.4(15) 

0.29 


17.6(15) 
0.29 



NOTE. — b : in unit of 10 13 erg cm 



Figure 4 compares the Suzaku XIS FI (XIS0+3) images in 
the 0.5 — lOkeV range, reconstructed by using the data during 
the first 20 ks (top) and the last 74 ks (bottom). The images 
are corrected for exposure and vignetting, and are non-X-ray 
(detector) background subtracted. The images were in addi- 
tion smoothed by a Gaussi an function with a = 0.' 17, follow- 
ing the procedure given in lMaeda et al.l (1201 lb . Note that the 
source is clearly detected i n the first 20 ks, but almost un- 
seen in the last 74 ks. Also ICheung et all (120121) argues that 
the 0.3-10 keV Swift XRT count rate of (6.1±1.8)xl0" 3 
counts s -1 is equivalent to a 0.5—8 keV flux of ~ 3.1 x 10~ 13 
ergs cm -2 s _1 , indicating a ~3x brighter source about one 
year prior to the Chandra observation, and providing further 
evidence of long-term X-ray variability in this source. 

Figure 5 summarizes the X-ray light curve of 1FGL 
Jl 31 1.7-3429 thus obtained during the AO-6 observation and 
compared with those obtained in the AO-4 observation. All 
the XIS0, 1, 3 count rates are summed in the various energy 
bands of 0.4-1 keV, 1-2 keV, 2-4 keV and 4-8 keV (from 
top to the bottom). Note that X-ray variability is clearly seen 
above 1 keV, but not in the 0.4—1 keV light curve. To see 



this more quantitatively, we calculated the nor malized excess 
varian ce a xxs (Vaughan et al. 2003; see also lKataoka et alJ 
120071) for the A04, A06 light curves and for the combined 
A04+A06 light curves. The cr xxs parameter is an estimator 
of the intrinsic source variance after subtracting the contribu- 
tion expected from measurement errors. As shown in Figure 
6 (top), <J 2 NXS is consistent with zero for the 0.4 — 1 keV 
band, whilst almost constant (within error bars) and signifi- 
cantly positive above 1 keVFI 

Figure 6 (bott om) presents the norma lized power spectrum 
density (NPSD; Havashi da et al.| [T998) which is a technique 
for calculating the PSD of unevenly sampled light curves. 
Note that data gaps are unavoidable for low Earth orbit X-ray 
satellites like Suzaku. Since the orbital period of Suzaku 
is ~ 5760 s, Earth occultations create artificial periodic gaps 
every 5760 s in the data, even when the source is continuously 

10 The fractional root mean square variability amplitude, -F V ar is often 
chosen in preference to cr'j slxs , although the two convey exactly the same 
information. Here, we prefer &% X S s ™P'y because variability is not sig- 
nificant for the 0.4 — 1 keV light curve and negative values of o~ xx g are 
possible. 
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FIG. 8. — Comparison of the average X-ray spectra of 1FGL J 13 1 1 . 7-3429 
observed in 2009 (A04; red) and 2011 (A06; blue) Suzaku observations. 
The X-ray data points represent the weighted mean of XIS0, 1 and 3. In addi- 
tion, a bow-tie shows the best-fit parameters of the correspon ding Chand ra 
source, CXOU J131 145.7-343030, observed in 2010 ICheung et aU2012D . 



monitored. To calculate the NPSD of our data sets, we made 
light curves of two different bin sizes of 256 s and 5760 s fo r 
each of the A04 and A06 light curves dKataoka et al]|2001h . 
The NPSD calculated for each light curves is well represented 
by a steep power-law with P(f) oc /- 2 -00±o.32 (aC)4) and 
P(/) oc /- 2 04±0-44 (A06 ) Note that me am pii tu de of the 
NPSD is larger for the A04 (red) data than in the A06 (blue) 
data. This is due to large flare observed in the first 20 ks of the 
A04 observation (Figure 4) and is consistent with the larger 
value of <J%xs compared to variability during A06 observa- 
tion (Figure 6 top). 

Figure 7 shows the folded X-ray light curves from the 
Suzaku A04 (top) and A06 (bottom) data sets. The XIS0,1 
and 3 data are all summed in 0.4—10 keV band. The phase 
zero is defined as MJD 56010.76808 as in the optical folded 
light curves shown in Figure 3 (bottom). The vertical axis 
represents the source counts divided by the average source 
intensity in the frame. Although there is a regular exposure 
gap for phases <f> = 0.3—0.4 due to Earth occultations, there 
appears to be a clear excess a t phases <f> = 0.4—0.6 in A04 
data (see, also iRomanil d2012l) ). Moreover, we found sim- 
ilar structure in the A06 data, but at different phases of (j) 
= 0.6—0.8 with much smaller amount of excess. To check 
that this excess is not due to artifacts caused by periodic 
gaps due to Earth occultation, we also made the folded light 
curves o f the nearby X-ray source in the same field of view 
(src B inlMaeda et alj ( rloTlh . or CXOU Jl 3 1 147.0-343205 in 
iCheung et all (120121) : see Figure 4), whose X-ray flux is com- 
parable to the average X-ray flux of 1FGL J131 1.7-3429. The 
results presented as dashed lines in Figure 7 show no clear 
excess for the nearby source. This suggests that the X-ray ex- 
cesses in the folded lighturves are most likely due to orbital 
motion of 1 FGL J 1 3 1 1 .7-3429 itself. 

Finally, Figure 8 compares the unfolded X-ray spectrum for 
1FGL J131 1.7-3429 averaged over the A04 (red) and A06 
(blue) observations. The best model fits for both observa- 
tions consist of power-law continua with photon indices, F 
~ 1.3 — 1.7. Significant curvature (or deficit of photons be- 
low 2 keV) was observed only in the 2009 data, which was 
tentatively modeled by an excess value of Galactic column 
density Nn\ N-r is however consistent with zero for the 201 1 
data. The model fitting results are summarized in Table 2. Al- 



though the nearby X-ray source (src B in | Maeda et ail d201 ll) . 
or CXOU J 1 3 1 147 .0-343205 in lCheung et alJ(l2012l) : see Fig- 
ure 4) is located only 1.6' apart to the south, contamination 
from this source and other nearby faint sources is estimated to 
be less than 5 % for the applied region of interest of 1 ' . 

4. DISCUSSION AND CONCLUSIONS 

During the first (A04; 2009) and second (A06; 2011) 
Suzaku observations, we detected significant variability in 
1FGL J131 1.7-3429 characterized by repeated flaring activ- 
ity, with a timescale of ~ 10 ks. The variability is only 
clearly seen above 1 keV. The NPSD is well characterized 
by P{f) oc /~ 2 , as is the case for the X-ray variability 
of va rious classes of AGN including Seyferts and blazars 
(e.g.. lHavashida et al.|[T998l: IVauehan et alj|2003l) . The latter 
class constitutes the majority of ferm/-LAT sources but the 
X-ray variability timescales of blazars are in general some - 
what longer, typically ~ 1 day (e.g.. Kat aoka et all 12001). 
Such "red-noise" PSD behavior is also observed in X-ray 
(Galactic black hole and neutron star) binary systems, but on 
much shorter time scales. For example, variability as short 
as ~ 1 — 10 ms has been observed fo r the famous Galac- 
tic black hole source Cyg X-l (e.g.. iMeekinseFall [1981: 
lHayashida eTaTlll998l) . 

Significant variability has also been observed in the optical 
(g', Re, r'), where it is rather a quasi-sinusoidal flu x modu- 
lation with a 1.56 hr period, as recently reported by Romani 
(120121) . Moreover, we also found that the modulation profile, 
including the amplitude of modulation and peak intensity, has 
changed largely among the six nights of observations. The ap- 
parent modulation of magnitude observed in both the IAO and 
LOT data is quite similar to those observed in 1FGL J2339.7- 
0531, which is character ized by a 4.63-hr orbital period in 
optica l and X-ray data ( Romani & Shawl 1201 lb iKong et all 
120121) . Note that 1FGL J2339.7-0531 is now suggested to 
be a "radio-quiet" gamma-ray emitting black-widow MSP 
with a ~0.1 Mq late-type companion star, viewed at incli- 
nation i ~ 57°. Moreover, this compact object companion 
in 1FGL J2339. 7-0531 is strongly heated, with T e g varying 
from ~ 6900 K (superior con junction) to < 3000 K at mini- 
mum (iRomani & Shawll201 ll) . 

In this context, the spectral energy distribution of 1FGL 
J131 1.7-3419 from radio to gamma-rays (Figure 9, top) may 
provide some hints on the nature of this mysterious source. 
1FGL J131 1.7-3419 has been the subject of both radio pul- 
sar counterpart searches and blind searches for gamma-ray 
pulsations, but to dat e no pulsed emission has been detected 
dRansom et alJ r20Tll) . and the situation is quite similar to 
the case for 1FGL J2339.7-0531. There are no NVSS radio 
sources catalogued within the 2FGL ellipse of 1FGL J 1 3 1 1 .7- 
3429 down to 1 .4 GH z flux completeness limit of ~2.5 mJy 
dCondon et al. 111998b . which is given as an arrow in Figure 9 
(top). From the 0.1 — 100 GeV gamma-ray flux listed in the 
2FGL catalog, F 1 = F .i_ioogcV — 6.2xl0~ n erg cm~ 2 
s -1 , we obtain F y /F x ^ 300, F^/F R > 1.7 xlO 6 for 1FGL 
Jl 3 11.7-3429, where Fx and Fr are the X-ray and radio 
fluxes measured in 2—10 keV and at 1.4 GHz, respectively. 
Note, these values are almost comparable to those measured 
in 1FGL J2339.7-0531, with F y ~ 3.0X10" 11 erg cm 2 s" 1 , 
F^IFx ^ 150 and F 1 IF R > 8.7 x 10 5 . Although its flat spec- 
trum X-ray continua and relatively high 7-ray-to-X-ray en- 
ergy flu x ratio (of > 100) i s typical of the sub-class of blazars 
(FSRO: lNoiariet^l2m2h . the non-detection in the radio as 
well as its curved gamma-ray spectrum seems to strongly dis- 



X-ray and optical observations of 1FGL J131 1.7-3429 
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FIG. 9.— (top) Broadband spectrum of 1FGL J131 1.703429. The X-ray 
data are as presented in Figure 8. The gamma-ray data points are taken from 
the 2FGL catalog (Nolan et al. 2 0121) . The radio upper limit of 2.5 mjy at 1.4 
GHz is taken from the NVSS catalog (Condon et al. 1998). The optical and 
UV data represents Swift UVOT data (see Maeda et al. 201 1 and Cheung et 
al. 2012). (bottom) Close-up of the broadband spectrum in the optical-UV 
range. Epochs 0—6 (1200 s each) correspond to time ranges defined in Figure 
2. Note the optical-UV spectrum is well fitted by a blackbody model of kT 
~ leV with a radius of the emission volume of ~ 1 .5 X 10 4 rfkpc km. 
where d^pc is distance to the source in unit of 1 kpc. 



favor an AGN association for 1FGL J131 1.7-3429. 

To further support the MSP association of 1FGL J1311.7- 
3429, Figure 9 {bottom) shows a close-up of the SED of 
1FGL J131 1.7-3429 in the optical/UV bands, reconstructed 
from IAO and Swift UVOT data. Ob served magnitudes wer e 
converted to flux density based on Fukug itaet al. I d 19951) . 
Red dashed line indicates a tentative fit with a blackbody 
model of T ~ 1.2 xlO 4 K (or kT ~ 1 eV) assuming an 
emission volume radius of R^ ~ 1.5 xlO 4 dkpc km, where 
dkpc is the distance to the source in units of 1 kpc. There- 
fore, the optical/UV spectrum of 1FGL J131 1.7-3429 seems 
compatible with what is expected from a companion star of a 
radio-quiet MSP like 1FGL J2339.7-0531. A slightly higher 
temperature than 1FGL J2339.7-0531 (kT ~ 0.3-0.6 keV; 
iRomani &~ Shaw 201 U) may indicate x2 smaller orbital ra- 
dius, ~ 5xl0 5 km, for the 1FGL J1311.7-3429 binary sys- 
tem, assuming that the mass of the companion star is 0. 1 Mq 
and the pulsar spin-down luminosity is L ~ 10 34 erg s _1 (pa- 
rameters suggested for 1FGL J2339. 7-0531 binary system; 
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FIG. 10. — (top) Comparison of the 2FGL Variability Index versus Cur- 
vature Significance for associated sources (AGN: blue, PSR: green) and 
unassociated sources (red), (bottom) Comparison of the 2FGL Photon In- 
dex versus 0.1 — 100 GeV Flux for the associated sources and unassociated 
sources. In both panels, a separation between the AGN and PSR populations 
is evident. Note that 1FGL 11311 .7-3429 is situated in the typical PSR region 
in the top panel, whilst at the boundary of AGN/PSR sources in the bottom 
panel. 



IRomani & Shawl [201 lb . A change of the modulation pro- 
file observed in the r'-band may be accounted for by rapid 
changes in the companion star temperature, but this remains 
uncertain. 

Since both the optical and X-ray light curves in the case of 
1FG L J2339.7-053 1 clearly exhibits a 4.63-hr orbital modula- 
tion (IRomani & Shawll201 it iKong et al.ll2012l) . the detection 
of periodicity in the X-ray light curve of 1FGL J131 1.7-3429 
is also likely. In fact, the folded X-ray light curve exhibits 
an exc ess feature around <j> = 0.5. As indicated by iRomanil 
(2012), this is presumably pulsar superior conjunction, al- 
though both the normalized intensity as well as the phase peak 
appear to have changed substantially between the 2009 (A04) 
and 201 1 (A06) observations. We therefore expect that X-ray 
variability consists of at least two different components - one 
associated with the binary motion as for the optical data, and 
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the other is rather random fluctuation well represented by the 
NPSD of P(f) oc f~ 2 , whose physical origin is still unknown 
but possibly related with perturbations associated with shock 
acceleration. 

Such flaring X-ray variability has not yet been observed 
for 1FGL J2339.7-3429, but solely in 1FGL J13 11 7-3429 . 
In a review of X-ray emission from MSPs, Zavlin (2007) 
describes three primary sources of X-ray emission: (1) 
intra-binary shock, (2) the neutron star (NS) itself, and (3) 
pulsar wind nebula o utside the binary system (see, also 
Archibald et alJ I2010I) . In fact, thermal emission of kT 
~0. 1 —0.2 ke V is often observed from MSPs, which is thought 
to ari se from the surface of the NS and is steady with time 
(e.g.. lMarelli et al.ll201 U iMaeda et al.ll201 ll) . More recently, 
thermal emission of kT ~ 0.1 keV was also detected from 
1FGL J2339.7-0531 (Kong et al. 2012, in prep), which is 
again steady despite large flux modulation associated with bi- 
nary motion being observed above 2 keV. Therefore, the fact 
that X-ray variability is not clearly seen below 1 ke V may sug- 
gest there could be some contribution from the surface of the 
NS, although from the spectral fitting, it is not statistically sig- 
nificant. Then the variable, hard X-ray emission could arise 
from the intra-binary shock rather than the nebula because 
variability as short as ~ 10 ks is unlikely to originate from 
the extended pulsar wind nebul aH Such a shock could re ad- 
ily produce gamma-ray emission dArons & Tavani 111993b . If 
localized, it could easily account for the orbital modulation 
as seen in the X-ray emission from 1FGL J2339.7-0531. But 
if material leaving the companion star, either through Roche- 
lobe overflow or a stellar wind is non-uniform or patchy, we 
might expect random flaring activity as we see in the X-ray 
data of 1FGLJ131 1.7-3429. 

We can also speculate on the nature of 1FGL J131 1.7-3429 
based solely on the gamma-ray pr operties, an approach al - 
ready applied to the 1FGL unlDs in lAckermann et aD (120121) . 
Figure 10 {top) presents a comparison of the 2FGL asso- 
ciated (either AGN (blue) or PSRs (green)) and unasso- 
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ciated sources (red) in the Variability.index and Sig- 
NIF.CURVE plane. Apparently, 1FGL J131 1.7-3429 is situ- 
ated in the typical PSR regions of this diagnostic plane. Sim- 
ilarly, Figure 10 (bottom) plots the distribution of PSRs and 
AGN in the PHOTON INDEX versus F .i_i oGeV plane. In 
this case, 1FGL Jl 31 1.7-3429 is at the boundary of AGN and 
PSR sources, so is still consistent with a PSR association. For 
comparison, we also plot the gamma-ray parameters for 1FGL 
J2339.7-0531. In both panels, the gamma-ray properties of 
1FGL J131 1.7-3429 and 1FGL J2339.7-0531 are quite sim- 
ilar. Again, this supports the idea that 1FGL J131 1.7-3429 
is a black-widow system and may be a second example of a 
"radio-quiet" MSP after 1FGL J2339.7-0531. 

Finally, if the rapid X-ray flaring variability observed with 
Suzaku may be due to inhomogeneity of shock material 
and/or rapid changes in the beaming factor, this could be ex- 
pected to occur also in gamma-rays, as suggested by a smooth 
connection of the spectrum between X-ray and gamma-ray 
energies. Moreover, we may see correlated variability also 
in the optical, which may be related to the change of modu- 
lation profile in the r'-band magnitudes we see in Figure 3. 
Unfortunately, such fast variability is difficult to observe in 
gamma-rays as we referred to one-month binned light curve 
in Section 1, despite the excellent sensitivity of Fermi-LPR. 
The low gamma-ray statistics also make it difficult to run a 
cross-correlation in order to measure any possible correlated 
optical, X-ray and gamma-ray variability. Further continuous 
investigation is necessary to confirm the origin of "variable" 
X-ray emission observed in the 1FGL J131 1.7-3429 system. 
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Variability has been observed in some pulsar wind nebula in X-rays 
and gamma-rays, but with typical timescales of a week to months (e.g., 



Pavlov et al. 2001; Abd o et aU20"TTl) . 



